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RESOLVED SHEAR STRESS INTENSITY COEFFICIENT 


AND FATIGUE CRACK GROWTH IN LARGE CRYSTALS 


INTRODUCTION 


Fatigue crack: propagation by cyclic plastic deformation on a slip 
plane early in the fatigue life is known as Stage I crack growth [1]. 

In Stage I crack growth in a polycrystal, the crack plane is a slip 
plane and it also coincides with the plane of the maximum shear stress. 

Neumann [2,3] and Vehoff and Neumann [4] have found that when a 
tensile stress was applied to a single crystal, a fatigue crack grew by 
the shear decohesion on two intersecting conjugate slip planes inclined 
to the tensile axis. The local small crack increments all took place 
on these two conjugate slip planes, while the macroscopic crack plane 
was perpendicular to the tensile axi3 . 

Even during stage II crack growth in a polycrystal, it has been 
generally agreed that a fatigue crack grows by the shear decohesion 
mechanism. Based on the shear decohesion process, Kuo and Liu [5], 

Yang and Liu [6], and Liu [7] have proposed an unzipping fatigue crack 
growth model. Their calculated fatigue crack growth rate agrees very 
well with the measured stage II growth rate and striation spacing. 

Extensive studies of fatigue crack growth have been conducted on 
Ni-base single crystals [8,9,10,11,12], and other aluminum-base single 
crystal alloys [13,14] as well as pure aluminum single crystals 



[15,16] . 


All of these studies have shown that fatigue crack growth is 


highly sensitive to the orientation of the crystal, and that the crack 
plane i3 crystallographic and follows a single slip plane or a 
combination of several slip planes. 

Since shear decohesion on a slip plane is caused by dislocation 
motion, the resolved shear stress acting on the active slip plane ahead 
of a crack tip must be resposible for the propagation of the fatigue 
crack as suggested by many researchers [10,13,15] . Several models have 
been proposed as the mechanisms of Stage I shear crack growth 
[9,15,17,18] . 

Recent work by Chan et al. [12] on fatigue crack growth in MAR-M200 
single crystals revealed that fatigue crack growth is mainly 
crystallographic cracking on {ill} planes. At room temperature, 
fatigue crack growth rates were found to be orientation dependent. 

Based on the elastic energy release rate associated with the mixed mode 
fracture, the effective stress intensity factor range AK eff was used to 
correlate crack growth, which incorporates K x , K IX , and K IIX for an 
anisotropic crystal. 

For fatigue crack growth in a single crystal or in a large grain, 
the anisotropy of plasticity must be taken into consideration. Hence, 
the stress intensity factors alone may not be sufficient for 
correlating with crack growth rate. In addition, the fracture plane of 
Stage I growth is crystallographic and often inclined not only to the 
loading axis but also to the specimen broad surface. Therefore, the 
mode of cracking is mixed, consisting of Mode I, II, and III 
components. A successful study must take both of these factors into 
account . 
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In our previous study of fatigue crack growth in coarse grain Al 
7029 aluminum alloy [19], the fracture surfaces were planar areas 
parallel to either {ill} or {100} planes. The {ill} crack surfaces 
were planar and shiny. They were formed primarily by shear decohsion 
on a single-slip-plane. The { 1 00} crack surfaces show "pine tree" 
morphology presumably formed by shear decohesions on two sets of 
intersecting slip planes. 

In this study, fatigue crack growth tests were carried out on large 
grain Al 7029 aluminum alloy. The crack path is usually very 
irregular. The finite element method is used to calculate the stress 
field near the tip of a zigzag crack. The resolved shear stresses on 
all twelve slip systems are computed, and the resolved shear stress 
intensity coefficient, RSSIC, is then defined. RSSIC is used to 
analyze the irregular c rack path, and it is used to correlate with the 
rate of single-slip-plane shear crack growth. 


EXPERIMENTAL PROCEDURES 


The material used in this study was Al 7029 aluminum alloy 
received as extruded 1/4" thick plates. The chemical composition in 
weight percent was 1.6%-Mg, 4.6%-Zn, 0.7%-Fe, and 0.05%- Si. The alloy 
was of high "purity" with very little impurities other than the specified 
chemical composition to promote grain growth. 

The plate was annealed at 540° C for about 50 hours followed by 
water quenching. The average grain size obtained after heat treatment 
was about 1 cm in diameter, with some exceptionally large grains 


3 



reaching 2.5 cm in diameter. 


Compact tension specimens were fabricated from the heat treated 
plates. The dimensions and a photograph of the specimen are shown in 
Figures 1 and 2. The outline of the grain boundaries can be clearly 
seen. 

Prior to the fatigue test, all specimens were mechanically 
polished and etched to reveal the grain boundaries. Fatigue crack 
growth tests were conducted under load-controlled conditions. The 
frequency was either 10 or 20 Hz at R - 0.1. The crack growth was 
monitored with a travelling microscope with a magnification of lOOx. 
The orientations of the crystals at the crack tips, as well as the 
fracture surfaces, were determined with X-ray Laue method. 

RESULTS AND DISCUSSIONS 

A. Orientations of the crack surfaces 

Figure 3a . i is a photograph of the broken specimen, R1 . Part of 
the fracture surface lies on a slip plane and part does not. The 
fracture area A is a non-slip-plane fracture surface, while area B is 
parallel to a {ill} slip plane. C designates the final area of 
fracture. The narrow strip bounded by dashed lines in the top sketch 
of Figure 3a . ii is the fracture surface. Fracture areas A, B, and C 
are also marked within the strip. Both of the side surfaces of the 
specimen are unfolded upward; the grain structure of the specimen is 
also shown in the figure. 

The lower sketch of Figure 3a. ii is the profile of the fracture 
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surface viewed from the broad side of the speciman on which the crack 
growth was monitored. Since area A is nearly perpendicular to the side 
face, area A becomes a line in the lower sketch. On the other hand, 
areas B and C are inclined toward the viewer, so both areas can be seen 
in the sketch. Similar sketches for other specimens are presented in 
Figures 3b through 3e. 

The macroscopic crack surfaces can be classified into two 
categories: slip-plane crack surface and non-slip-plane crack surface. 

A crack surface often changes its orientation from one grain to the 
next, and a crack may switch from slip-plane cracking to non-slip-plane 
cracking and vice versa, forming a zigzag crack path. The orientations 
of these crack planes were determined by the X-ray Laue photographs 
taken either from the crack surface or from the side face of the 
specimen. Slip-plane cracking can be easily identified from their 
bright appearance and the straight and sharp intersections with 
specimen sidefaces. 

Slip-plane cracking is caused by shear decohesion on a single slip 
plane, while non-slip-plane cracking is a result of alternate shear 
decohesion on two or more different slip planes. Shear decohesion is 
due to dislocation motion, which is controlled by the forces on the 
dislocations at a crack tip. The force on a dislocation is directly 
related to the resolved shear stress on the slip plane. Therefore, 
the rate of shear decohesion must be related to the resolved shear 
stress "intensity” at a crack tip. 

The c rack tip stress field will be calculated by the finite element 
method for the zigzag crack path of each specimen. With the 
orientation of the crystal at a crack tip determined by the Laue X-ray 
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pattern, the resolved shear stress field on each of the twelve slip 
systems at a crack tip is computed and its intensity evaluated. 

The active shear decohesion plane or planes must be the slip 
plane (s) with a high resolved shear stress intensity. If the resolved 
shear stress intensity on a plane is much higher than all the rest, the 
plane must be the primary plane of shear decohesion and the slip plane 
becomes the crack plane. On the other hand, if the resolved shear 
stresses on two slip systems are comparable, shear decohesion will take 
place on both of these two slip planes, and the macro crack plane will 
not follow a single slip plane. The orientations of the crack surfaces 
will be analyzed with the computed resolved shear stress intensities. 

If a crack surface is a slip plane, it is logical that the crack 
growth rate on the single slip plane will correlate with its resolved 
shear stress intensity. 

B. Finite Element Modelling of the Crack Tip Resolved Shear Stress 

The sketches in Figures 3a to 3e indicate that crack paths are very 
irregular and zigzag. The elastic crack tip stresses of such zigzag 
crack are calculated by FEM. A two-dimensional plane model of a 
through crack is used to simulate the actual crack. Both Mode I and 
Mode II crack tip fields are calculated. 

The real crack surface is three-dimensional and often inclined to 
the broad surface of a specimen. The inclination of the crack surface 
induces Mode III crack tip field. The Mode III crack tip field is 
obtained by an approximation. 

Once the Mode I, II, and III crack tip fields are obtained, the 
resolved shear stresses on various slip systems of the crystal at the 
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crack tip can be calculated. The resolved shear stress field of a slip 
system is defined by its intensity coefficient. 

Figure 4 shows the sketch of a crack in a compact tension specimen 
for the FEM simulation. Although a real crack may be zigzag, the c rack 
in the simulation consists of a machined notch and an inclined crack. 
The inclined crack is a two-dimensional through crack and the crack 
surface is perpendicular to the broad specimen surface. The 
inclination angle a is defined as shown in the figure. 

The local coordinate with the origin at the crack tip is shown in 
the figure. The x'-axis is defined as the line along the inclined 
crack; and y'-axis, perpendicular to the crack line. The r and 0 are 
local polar coordinates, defined in the conventional manner. 

Figure 5 shows a typical finite element mesh for a compact tension 
specimen. The ABAQUS FEM program was used. The crack tip elements 
consist of eight layers. The plane eight-noded quadrilateral 
quadratic isoparametric element was used. The crack tip triangular 
elements were formed by collapsing one side of the quadrilateral 
element. The mid-node of a crack tip element was moved to the quarter 
point to take account of the elastic singularity of the crack tip 
field. 

To insure proper finite element modelling, the crack tip stress 
field of the conventional compact tension specimen with a straight 
crack was calculated. The calculated stresses at 0=0° were plotted as 
a function of r, the distance from the crack tip, and the results were 
compared with known solution. As shown in Figure 6, the finite element 
calculation yields very satisfactory results. 

In the case of an inclined crack, the crack tip stresses are 
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given by: (referring to Figure 4 for the definition of the local 

coordinates) 


’ i,j ' ' vib (9) + ai ' j, ' n (9) 


(i) 


Along the crack line, i.e. 9“0°, or along the x'-axis. 


°x'x’ 


<*y»y ' 


K I 

V2nr 


X x ' y ' 


Kn 

V2Jlr 


( 2 ) 

(3) 


Thus, the stress intensity factors, K x and K n can be evaluated from 


K I * <Ty'y« ^2jtr 


(4) 


K ri = ^x'y ' ^2Jtr 


(5) 


where Oyiyi and x x i y < are the computed local stresses along the crack 
line, which are obtained by transforming the calculated stress 
components of a inclined crack: into the x'-axes. 

Figure 7 shows an example of the log-log plot of Oyiyi and T x iyi 
v.s the distance from the crack tip, r. Figure 8 is a plot of K x and 
K n computed from Equations 4 and 5 as functions of log r. Notice that 
K l and K n are nearly constant at the distance very close to the crack 
tip. These asymptotic values of K T and K n are, therefore, taken as the 
norminal K r and K n for the crack geometry. 

The crack tip field for a through crack with a crack surface 
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inclined to the broad surface of a specimen is characterized by K T K n 
and K m . A two-dimensional FEM calculation can only give K x and K n . 
As a consequence, a procedure is adopted to obtain the approximate 
values of K z f K n and K m for the inclined crack. The detailed 
derivations are given in the Appendix. 

The approximation for the inclined crack was proposed by Pook [20] 
and Chan and Cruse [21] . The calculated Mode I and Mode II crack tip 
fields are projected onto the inclined crack surface. The equations 
for the K's in the case of three-dimensional inclined crack are given 
by 


K\ - Kj cos 2 P ( 6) 

K n = K n CO3 0 (7) 

K' ni * Kj cos P sin P (8) 


where K T and K n are stress intensity factors obtained from the two- 
dimensional finite element calculation; p is the angle of the 
inclination of the crack surface, defined as the angle between the 
normal of the crack plane and the y'-axis. 

For a three-dimension inclined crack geometry, the crack tip stress 
components are given by 




V2 7Cr 


CT i ’ ) ' , r ( 9 ) + 


^2kc 


CT i ' 1 ' , 11 ( ®) 


K', 


V2H 


111 ~ , 

CT i 'j ' ,m ( 0 ) 


(9) 


r and 0 are the local polar coordinates at the crack tip as shown in 
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Figure 4 . 


In an earlier study by Duquette et al. [10] on Stage I fatigue 
cracking in a Nickel-base superalloy single crystal, the glide forces 
on a dislocation lying in all possible slip planes were calculated. 

The calculation showed that crack growth is directly related to slip 
activities on slip systems having the highest resolved shear stresses. 

Similar procedures were followed by Nageswararao et al. [13] on 
Stage I fatigue propagation study in an Al-Zn-Mg alloy. Relative shear 
stress values on various slip systems at a given distance from the 
crack tip were computed. Fatigue cracks propagated preferentially in 
the slip system that experiences the maximum resolved shear stress. 

The resolved shear stress is given by [22] 


1 

t rss * k ^ij n j 


( 10 ) 


where and b are the Burgers vector and its magnitude; nj the unit 
normal vector of the slip plane; and the crack tip stress tensor 

field. Substituting Equation 9 into Equation 10, the resolved 
shear stress is 


Trss “ vlfc [ bl " ]I K i' k i'i ' K m ' f (6)] [ «=i - 3 •! T | nj»] (11) 

As defined in the Appendix, x^" are the principal axes of the crystal. 
C^Hji are the direction cosines between x^ M and x j ' . b^». is the Burgers 
vector; n j »» the unit normal vector of the active slip plane referring 
to the x^” coordinate system. The matrix containing K’s and f(0) is 
the stress matrix, whose component is defined by Equation 9. 

10 


l 



Equation 10 indicates that t RSS preserves the 1/Vr singularity, and 
the intensity of X RSS is dependent on the crystal orientation relative 
to the crack surface. For a given crystal orientation and crack 
geometry, the angle 0 is equal to 0 3 , the angle between the 
intersection of the slip plane with the broad surface of the specimen 
and the x'-axis as defined in Figure 4. x RSS is a function of r, 0 3 
and the values of K's. 

The intensity of X RSS is linearly proportional to the quantity 
RSSIC, which is defined as following 


RSSIC - [ b l «] [ C iWj •][ Kj 


'n ' 


^TTT • 


f (®s^] [ " j 1 ] [ n j ”] 


( 12 ) 


We call^the quantity the resolved shear stress intensity coefficient 
(RSSIC) . For a given slip system, X RSS is 


x 


R S S 


RSSIC 

V27Cr 


(13) 


The intensity of x RSS of this slip plane depends on the RSSIC 
value, which is analogous to stress intensity factor K. However, K is 
applicable to any angle 0, while RSSIC is applicable only to a specific 
slip plane. For different slip systems the values of RSSIC f s are 
different. RSSIC can be used to characterize the resolved shear stress 
of a given slip system and the forces on the dislocations of the slip 
system. Therefore, it can be used to characterize the slip behaviors 
of the slip system and the shear decohesion process at a crack tip. 

Two advantages in using RSSIC are: (1) the dependence of x RSS on r 

is eliminated; (2) the angle 0 S has a definite physical meaning, which 
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is directly related to the orientation of the slip plane. 

jn the following section, the detailed analyses are presented for 
the relationship between the resolved shear stress intensity 
coefficients and the fracture planes of the tested specimens. 

C. RSSIC and the Active Slip Plane 

It is well established that ductile fatigue crack growth is caused 
by shear decohesion [2,3,4,23,24,25]. Shear decohesion is caused by 
dislocation motion. Therefore, it is expected that the orientation of 
the crack plane must be related to the active slip plane (s) . A 
detailed quantitative study on the RSSIC and its relation to the active 
slip plane (s) will be very useful to advance the 3tudy of the shear 
decohesion process . 

Four speicmens were tested and analyzed. Crack surfaces of two 
specimens follow the {ill} slip planes, and another two specimens have 
non-slip-plane crack surfaces. 

1. Specimens with slip-plane fracture surface 

Specimens R2 and R3 are good examples of slip-plane fracture 
surface. The sketches of the crack surfaces are shown in Figures 3b 
and 3c. Figure 9 shows a photograph of the fractured specimens R2 and 
R3 . The slip-plane fracture surfaces of these two specimens are very 
shiny, smooth, and reflective to naked eyes. 

The orientations of the crystals at the crack tips were determined 
by x— ray Laue method. They are very close to each other. A 
gtereographic projection of the crystal is presented in Figure 10. The 
loading axis (LA) and the crack propagation direction (CPD) are marked 
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in the figure. The macroscopic crack plane of each of these two 
spcimens followed a single {ill} plane. 

For Specimen R3, the crack was initiated at the center of the 
machined Chevron notch because of the high tensile stresses there. 

After the crack initiation, it was observed that on one side of the 
specimen surface, the crack propagated vertically downward at first. 
Then the crack turned forward. On the other side of the specimen, 
however, the crack propagated straight forward, but slightly off the 
specimen center line. After the initial transition stage, the crack 
front shifted to the (111) slip plane, which is inclined to the 
thickness direction of the specimen. From thereon, the crack continued 
to propagate on this plane until final fracture. A similar behavior 
was observed on Specimen R2 . 

The stress state in the initial stage is too complicated 
to investigate. The analysis is made mainly for the major portion of 
the crack, i.e. along the (111) crack surface. 

First, we have to calculate the RSSIC of each slip system and to 
locate the active slip plane (s). Both specimens have the same crack 
angle (i.e. a-3°) and same crystal orientation with respect to the 
local coordinates at a crack tip. The relative orientation is given by 
the matrix of the direction cosines [Cj^wji]. With a, and the 

crack tip field, i.e. K r f , K r * T , and , known, the resolved shear 
stress intensity coefficients, RSSIC, for all 12 slip systems are 
calculated with Equation 12. The results of the calculation are 
tabulated in the first column of Table 1 for the projected crack length 
of a ' /W - 0.423, where a 1 * t + a cos a , as indicated in Figure 4. 

The solid line segment at the left hand side of the sketch at the 
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bottom of the first column in Table 1 is the notch. The next solid line 


segment to the right is the crack, which follows (111) plane with an 
inclination angle a * 3°. The dashed line segments indicate the 
orientations of other {ill} planes. 

From the table, one can see that the value of RSSIC on (111) [101] 
slip system (underlined) is much higher than all others. Therefore, 
this must be the dominant crack tip slip system. The orientation of 
(111) slip plane coincides with that of the crack surface. 

We also explore the possibility that the carck may turn to one of 
the other three slip planes. Additional calculations with (111), (111) 

and (111) crack increments are made. 

The values of RSSIC for the crack increment along the (111) slip 
plane are listed in the second column. The highest value takes place 
on the (111) [Oil] slip system. Therefore, once a crack follows the 
(111) plane, the crack will switch away from the (111) plane to the 
(111) plane. 

The calculations for the crack increment along the (111) -plane are 
listed in the third column. The (111) [Oil] has the highest RSSIC. 
However, the value is 3till much lower than that of the (111) [101] slip 
system. The RSSIC controls the dislocation motion. Therefore, the 
shear decohesion process takes place on the (111) -plane as observed. 

If localized softening takes place in the active slip band, a further 
concentration of the slip motion in this band will be promoted as 
suggested by Nageswararao et al. [8]. 

Figure 11 is an optical photograph of the fracture surface. Many 
parallel slip markings across the specimen thickness are clearly 
visible. The markings spread over the entire fracture surface until 
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the final area of failure. From the x-ray Laue photograph, these slip 


markings were identified to be parallel to the [101] direction, which 
is the intersection between the (111) crack plane and the (111) slip 
plane. It is interesting to observe that with the (111) crack plane 
(Column one), the RSSIC's of the (111) [101] and (111) [110] slip systems 
are also very high. These markings indicate that one or both of these 
two slip systems were also active. However, the shear decohesion on 
this slip plane did not make any significant contribution to the crack 
growth rate . 

The values of RSSIC for the crack increment along the (111) -plane, 
the fourth column, are much lower than those for the crack increments 
along the other three slip planes. Therefore, the (111) -plane will not 
participate actively in the shear decohesion process. 

Since Specimen R2 has the same crystal orientation and fracture 
surface as R3, the above calculations, analyses and conclusions are 
equally applicable. 

For a single active slip plane, the maximum resolved shear stress 
provides the primary driving force for dislocation motion, the shear 
decohesion process, and, thereby, the crack growth process. The 
correlation of the crack growth rate on such plane with the calculated 
range of RSSIC, ARSSIC, will be shown in a later section. 

2. Specimens with non-slip-plane fracture surface 

The immediate extension of the above observations and analyses is 
that if the values of the RSSIC on two or more slip systems are equal 
or comparable, these slip systems will be activated and shear 
decohesion can take place on these slip planes . 
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The model for ductile fracture and ductile fatigue crack growth 
proposed by Orowan [26], Neumann [3] r McEvily [15], Laird [23], Pelloux 
[24], and Tomkins [25] involves the alternate shear decohesion process 
on two intersecting conjugate slip planes passing through the crack 
tip. 

When the alternate shear decohesion process takes place on two 
intersecting conjugate {ill} slip planes in a FCC crystal, the 
macroscopic fracture plane will coincide with a {100} plane as observed 
by Neumann [3] in copper single crystals. The copper crystal was 
oriented so that the notch plane was parallel to the {100} plane, and 
two {ill} slip planes were symmetrical to the notch plane, i.e. one 
above and the other below the notch plane. See Figure 12. The shear 
decohesion process alternated between these two {ill} slip planes. 

While the microscopic crack plane always follows the {ill} planes, the 
resultant macroscopic fracture surface is parallel to the {100} plane, 
i.e. a non-slip-plane fracture surface for the copper crystal. 

The distinction between the slip-plane and the non-slip-plane 
fracture surfaces can be readily seen from the morphologies of the 
fracture surfaces of Specimens R4 and R5 in Figures 13 and 14, 
respectively. The common features of the fracture surfaces of these 
two specimens are very rough surface, divergent river lines originating 
from the notch tip, and the lack of macroscopically well-defined slip 
markings . 

The stereographic projections of these two cystals are shown in 
Figures 15 and 16. The macroscopic crack planes were also determined 
by the Laue X-ray method. These planes are very close to { 1 00 } -planes . 
The stereographic projection of the crack plane of Specimen R5 is shown 
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in Figure 17. Notice that the center of the projection is very close 
to the {100} pole. Both crack planes were nearly perpendicular to the 
specimen broad surfaces, i.e. P - 0°. The angles of inclination, a, 
are 10 degrees and 33 degrees for R4 and R5, respectively. 

For Specimen R4, four sets of calculations are made. The first set 
is for the basic inclined crack at a - 10°. Three additional sets of 
calculations are made for the crack configurations shown at the bottom 
of each column. Each crack configuration has a segment as the notch 
and a segment as the inclined crack at a - 10° followed by additional 
segment (s) as shown in the sketches. 

The first column is the calculation of the crack tip stress field 
for the projected crack length of a'/W - 0.374. The initial solid line 
segment at the left hand side of the sketch at the bottom of the first 
column in Table 2 is the notch. The next solid line segment to the 
right is the crack, which has an inclination angle a - 10°. 

After the crack tip stress field is calculated, the values of RSSIC 
of the 12 slip systems are then computed. The results are listed in 
the first column of Table 2. In the table, the higher values of RSSIC 
are underlined. 

The slip sy 3 tem 3 (111) [101] and (111) [110] experience much higher 
resolved shear stresses. The (111) plane is located above the crack 
plane; while the (111) is below the crack plane. This suggests that a 
crack may grow alternately on these two slip planes and the resultant 
macroscopic fracture plane is somewhere in between these two planes. 

The actual crack plane nearly bisects the angle between these two 
active slip planes. 

Since the RSSIC of the (111) [110] slip system is nearly 20% higher 
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than that of (111) [101] slip system, it cannot be ruled out that the 
crack may propagate on the single slip plane of (111) . 

Assume that the crack will grow on the (111) slip plane. A finite 
element calculation is made for the new crack geometry. After the 
initial crack segment at a * 10°, a small crack increment along the 
(111) -plane is added as shown by the sketch at the bottom of the second 
column . 

The second column in Table 2 shows the calculated RSSIC of all slip 
systems based on this new crack geometry. Ralatively high RSSIC values 
are retained on these two slip systems. This implies that there is no 
substantial difference in the driving forces for dislocation motion 
between these two slip systems. Even if the crack grows on the (111)- 
plane, the crack may switch back to the (111) -plane. Therefore, the 
initial conclusion that shear decohesion will take place on (111) and 
(111) -planes is still valid. 

Additional calculation i3 made by assuming that the crack further 
switches from the (111) -plane to (111) -plane. An additional crack 
segment along the (111) -plane is added as shown by the sketch at the 
bottom of the third column. The calculated RSSIC value is listed in 
the third column. Again, the RSSIC values on these two slip planes 
remain higher than those of the others. 

It should be pointed out that the sequence in which the crack 
switches between (111) and (111) planes does not affect our conclusion. 
This is depicted by the calculations by assuming that the crack will 
switch from the inclined crack at a = 10° onto the (111) -plane 
directly. The result is listed in the fourth column of the table. 

These calculations demonstrate that no matter on which slip plane 
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the crack propagates, the RSSIC on the (111) and (111) slip plane3 have 
the highest values. The (111) plane is above the crack plane and the 
( 111 ) is below the crack plane. Therefore, the shear decohesion 
processs will alternate between two planes, and results in a non-slip- 
plane fracture surface very close to the (100) -plane. 

As determined from the stereographic projection of the crystal, the 
(111) and (111) planes intersect the crack plane on the broad surface 
of the specimen by -54 degrees (0 3 — 54°) and 59 degrees (0 3 «59°) , 
respectively. The crack plane deviates from the bisectional plane of 
these two 3 lip planes by 2.5°, i.e. ~ (59°-54°) . 

Poj- specimen R5, the inclined crack angle, ct, is 33®. Four sets of 
calculations were made for the crack geometries shown by the sketches 
in Table 3. The calculations were made for the initial projected crack 
length of a'/W — 0.32. 

First column of the table shows the values of RSSIC for the crack 
with an inclined angle of 33°. The two highest values are RSSIC - 2.01 
for the ( 111 ) [ 110 ] slip system and 1.64 for ( 111 ) [ 110 ] slip system. 
Additional calculations were made for the crack increments along the 
(111), (ill) and (111) -planes as shown by the sketches below. 

The higher values of RSSIC are underlined. The results indicate 
that if a crack propagates along (111) -plane, it will continue to grow 
on this plane. Therefore, it becomes and slip-plane cracking. 

However, it also indicates that once a crack switches to either (111) 
or (ill), the RSSIC values for the (111) -plane slip systems are much 
lower. Therefore, once switches to the (111) and (111) -plane slip 
systems, the crack will not be able to go back to the (111) -plane. 
There are a number of dislocation barriers in a crystal to cause cross 
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slip from (111) to (111) or (111) . (Ill) is above the crack plane and 
(111) is below. The crack will grow by the shear decohesion process 
alternately between (111) and ( 111) -planes . The crack plane bisects 
these two slip planes and it is close to the (010) plane. 

D. Growth Rate of Single Slip-plane Shear Crack 

The discussion thus far have clearly indicated that the resolved 
shear stress coefficient is capable of characterizing slip behaviors on 
specimens with either slip-plane fracture surface or non-slip-plane 
fracture surface. Since shear crack growth is a shear induced slip 
process, the parameter directly related to localized shear stress 
acting on specific slip plane, RSSIC, is a "natural” candidate as the 
controlling factor. In terms of dislocation mechanics, the glide force 
is the resolved shear stress acting on the slip plane in a slip 
direction, which provides a driving force for dislocation motion. 
Therefore, the maximum resolved shear stress intensity coefficient, 
RSSIC, on the crack plane (also a slip plane) is proposed as a 
parameter to characterize the shear crack growth. 

Figure 18 shows the maximum RSSIC on the slip system (111) [101] of 
Specimen R3 v.s. crack lengths, a'/W. The crack growth rate, da/dN, 
is plotted as a function of ARSSIC in Figure 19. 


DISCUSSIONS 


The finite element calculations for the crack tip stress fields are 
made for an isotropic solid. However, the aluminum crystal is 
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anisotropic with a cubic symmetry and the inclined crack needs a three- 
dimensional model. However, the calculations made in this study should 
be able to include the essential features of the crack tip stress field 
and to cover the primary effects on the crack tip slip systems. 

A crack tip slip system will be activated and shear decohesion will 
take place only if its resolved shear stress intensity exceeds a 
certain critical value. The resolved shear stress intensity is 
linearly proportional to the applied nominal K-value. For a specific 
crack geometric configuration, at a low level of the nominal K, there 
might be only one slip system activated, while at a higher level of the 
nominal K, two or more slip systems might be activated. The detailed 
shear decohesion process and the crack surface morphology depends on 
the orientations of the active slip systems and their relative 
contributions to the overall crack growth process. 

It is reasonable to expect that shear crack growth should correlate 
well with RSSIC. The crack growth rate correlation with the calculated 
RSSIC, Figure 19, is primarily for Mode III shear crack with a minor 
Mode II component . It remains to be seen whether the same correlation 
is equally applicable to Mode II, Mode III, and all of the combined 
shear mode crack growth. 

Fatigue crack growth on two or more slip planes might be related to 
a combination of the RSSIC* s of the active slip systems. It is evident 
that a large amount of information is needed to answer all of these 
questions. It is hoped that this promissing preliminary study will 
lead to a large effort in the basic study of the fatigue crack growth 
mechanism. 
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SUMMARY 


1. Fatigue crack growth is caused primarily by shear decohesion at a 
crack tip, which is the result of the dislocation motion and is 
controlled by the resolved shear stresses on the active crack-tip 
slip systems. 

2. Fatigue tests were conducted on A1 7029 large grain poly-crystals. 
The crack paths were generally irregular and zigzag. The crack tip 
stress field is calculated with the finite element method. The 
crack tip resolved shear stresses are computed, and the resolved 
shear stress intensity coefficient, RSSIC, is defined and the 
values for each test specimens are evaluated. 

3. When the RSSIC on a single slip system is much higher than all the 
others, the crack will follow a single slip-plane. 

4. When the RSSIC* s on two conjugate slip systems are comparable, a 
crack will grow in a zigzag manner on these planes, and the 
macro-crack-plane bisects these two active slip-planes. 

5. The maximum RSSIC on the most active slip system is proposed as 

a parameter to correlate with the shear fatigue crack growth rate 
in large crystals. 
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APPENDIX 


A compact tension specimen with its coordinate systems is shown 
schematically in the figure below. Xj are the coordinates with the 
origin at the notch tip and x L along the notch line. The origin of the 
Xj' coordindates is at the crack tip. In the two-dimenional plane, x’ 
is along the crack line and x’ is normal to the crack line. 



Since the crack plane is inclined to the loading axis and it is 
also inclined to the x'-xj plane, the crack tip stress field is a 
combination of Mode I, II, and III. Thu3, the crack tip stresses, 
referring to the x t ' -coordinates, are prescribed by Kj , , and K^ n . 

With a through crack and with the crack plane perpendicular to the 
x * -x 2 ‘ plane, the two-dimensional K T and K n can be calculated with FEM. 
However, if the crack plane is inclined to the x' -x' plane, in 
addition, Kj U exists also. By resolving the applied stress onto the 
crack plane, Pook [20] and Chan and Cruse [21] approximated the crack 
tip field for the inclined crack from the two dimensional crack tip 
field of Kj and K n : 


K j = Kj cos 2 P A1 

K ' n = K n cos P A2 
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K ' m ■* K r cos p sin p 


A3 


P is the angle of the inclination defined in the above figure. 

The crack tip stresses in the coordinates system are therefore 
given by [27] 

i 

k t 
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<T. ». • - . :■ cos: 
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Now, let x” be the local coordinates at the crack tip, and x[ f 
coincide with the principal axes of the crystal, i.e. x” in the 
direction of [100] ; x’\ [010] ; and x", [001] . At any local point, the 
stress components in the direction of the principal axes of the crystal 
can be denoted as O i n^n They are given simply by transforming <J t i ^ « 
into the directions of the crystal axes: 


"i 


[ C l ft j <] [ <?! • j '] [ C i " j '] 


A5 


where [C l .*-,»] is an 3x3 matrix, whose componenets are the direction 
cosines between x.J and x” axes. [<T 1 •] is an stress matrix, whose 
components are given by Equation A4 . 

The resolved shear stress is given by [22] 


TT b i CT u 


A6 


where b 1 is the Burgers vector, and is the unit normal vector of the 
pertinent slip plane. For a FCC crystal, b 1 is along <110>, and n ^ is 


24 


I 



along <111>. Substituting Equation A5 into Equation A6, the resolved 
shear stress in matrix form is 

^ R S S ~ ’*] [ w j *] [ 1 j 1 ] [ ^1 ” j 1 ] [ n j ”] ^ 

Since the components of [ O i i ^ t] are a function of K x , and 0, 

at a given local point near the crack tip, [ a A i 3 i] can be abbreviated 

as 

[ o^.] - -j==- [ Kj t k^, K; n , f (0)] . as 

\27Cr 

Hence, the resolved shear stress can be written as 

X R s s " t [ ”] [ ” j 1 ] [ K i * K n 9 9 ^ ” j 1 1 [ n j "] • 

\27Cr 

b t , n^ and 0 are related. Once b L and n^ are given, 0 is determined. 
We define the resolved shear stress intensity coefficient (RSSIC) as 

RSSIC - [ b^H C lWj »][ Kj , K£, K^, f (0)][ C lWj .] T [ n 3 ..] , A9 

such that 


RSSIC 

Trss “ ^ ‘ 
RSSIC is in the unit of MPaVm or KSlVin . 


A10 
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Table 1. RSSIC Values of 12 Slip Systems for Specimen R3 


Slip System 

(111) [10T] 
(111) [Oil] 
(111) [lTO] 

(Hi) [101] 
(Hi) [110] 
(Hi) [Oil] 

(Til) [OlT] 
(Til) [110] 
(Til) [101] 

(llT) [lTO] 
(llT) [101] 
(llT) [Oil] 


Crack Geometry 



RSSIC : 

MPaVm 


12SZ 

0.913 

3.342 

1.477 

3.599 

1.201 

1.177 

0.802 

1.667 

0.289 

2.165 

2.278 

4.357 

1.861 

1.902 

MU 

3.875 

0.913 

0.909 

2.215 

0.482 

0.949 

2.811 

1.599 

0.741 

4,273 

4.487 

1.413 

1.845 

3.677 

0.553 

1.371 

2.586 

0.597 

3.934 

2.784 

0.128 

1.783 

0.380 

3.216 

3.551 

1.614 

3.358 

2.273 

3.679 

3.397 

3.737 

0.943 

crack (1*^) 

(Hd'Vh) 

, <UT) 
cra< *j^mn 
(lTl)i'' (Til) 

— T\ 

(lTl)*'(Tii) 

SiniD 
_ t'-di 

( iTi)N(Tu: 
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Table 2. RSSIC Values of 12 Slip Systems for Specimen R4 


Slip System 


RSSIC 

: MPaVm 


(111) [101] 

0.460 

0.682 

0.120 

0.102 

(HI) [OlT] 

0.358 

0.007 

0.662 

0.566 

(111) [HO] 

0.818 

0.689 

0.542 

0.464 

(Hi) noli 

m 

2.224 

2.746 

2151 

(Hi) [110] 

1.348 

1.940 

0.640 

0.545 

(lit) [Oil] 

0.949 

0.284 

3.387 

2.898 

(Til) [Oil] 

0.287 

1.479 

2.050 

1.753 

(Til) [110] 

0.829 

1.852 

0.043 

0.038 

(111) [101] 

1.116 

0.372 

2.093 

1.791 

(llT) [lTO] 

211A 

2J31 

3,044 

2M& 

(ill) not] 

1.474 

0.946 

1.289 

1.104 

(HI) [Oil] 

1.340 

1.352 

1.755 

1.502 


Crack Geometry 



(tii) 

--(Til) 

'(iTt) 



(HT) (in) (111) (in) 

c rack J^'( Tll) crack/\/'^ Tu > 

\(lTl) \oT„ 
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Tabic 3. RSSIC Values of 12 Slip Systems for Specimen R5 


Slip System 


RSSIC 

: MPaVm 


(HI) HOT] 

0.532 

0.392 

0.568 

0.962 

(111) (Oil] 

1.177 

0.994 

1.556 

0.750 

(111) [iTO] 

0.645 

1.386 

0.987 

0.212 

(111) [101] 

1.116 

0.830 

2,532 

2,682 

(111) [HO] 

1.641 

0.798 

0.574 

2.084 

(Hi) [Oil] 

0.526 

1.628 

2.008 

0.598 

(Til) [Oil] 

1.295 

2.261 

1.191 

0.677 

(Til) [110] 

2M1 

1.446 

0.467 

0.167 

(Til) [101] 

0.716 

0.815 

1.659 

0.510 

(111) [ITO] 

0.963 

1.740 

0.155 

0.921 

(llT) [101] 

0.235 

0.468 

1.939 

2.285 

(11T) [Oil] 

1.198 

1.273 

1.784 

1.363 


Crack Geometry 


(Hi) 

1/ (no 

c rackyK: (nl ) 

’(til) 


|/ m > 
1/ (HD 


crack y 


on) 

x;-(lll) 
<1 IT) 
(III) 


crack 

— S (III) 


I) 



(III) 


l) 

(ill) 
(1 IT) 
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Fig. 1 Dimensions of compact tension specimen. 
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Fig. 2 The photograph of a compact tension specimen with large grains. 
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(3a.i) Photograph of the fractured specimen Rl. 


Fig. 3 Photograph and sketches of fracture surfaces and grain structures of specimens. 
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(3a.ii) The sketch of the fracture surface of Specimen Rl. 
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(3c) The sketch of the fracture surface of Specimen R3. 
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(3e) The sketch of the fracture surface of Specimen R5. 
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Fig. 4 A compact tension specimen with an inclined crack 
and its coordinate systems. 
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Fig. 6 The comparison of the FEM calculations of the 





(ly-Ddw) 


i 


44 


Fig. 8 The asymptotic values of the Kj and Kjj as r approaches to the crack tip. 



Fig. 9 Photograph of Specimen R2 and R3 with (111) slip-plane fracture surfaces. 
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Fig. 14 The optical photograph of the fracture surface of R5. 
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Fig. 17 The stereographic projection of the crack surface of Specimen R5. 
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Fig. 18 RSSIC as a function of the crack length for the highest RSSIC slip system, (1 1 1)[ 101], of Specimen R3. 
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